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OUTLINE 


Computational Geometry - how it fits in 
Survey - recent work 

A Computational Geometry Approach - current work 


COMPUTATIONAL GEOMETRY 


The design and analysis of algorithms and data 
structures for the solution of geometric problems. 



WHY COMPUTATIONAL GEOMETRY 


Complexity 

Bounds 

Robustness 


"This program takes 2 minutes to generate a 
grid for model X on workstation Y." 


Questions: 

Does the program always generate a grid? 

How does the number of grid cells affect execution time? 
What can be said about grid quality? 


293 


"CT-Notation 


A function T(n) is 0(f(n)) is there exist constants 
c and no such that for ail n>n 0 , T(n)<c f(n) 

Delaunay Triangulation - 0(n log n) 

Shamos and Hoey - Divide and conquer 
Fortune - Sweepline 

Guibas, Knuth, Sharir - Randomized incremental 


OPTIMALITY CRITERIA 

The Constrained Delaunay Triangulation 
minimizes the largest circumcircle 
minimizes the largest min-containment circle 
maximizes minimum angle 

lexicographicaly maximizes list of angles, smallest to largest 
minimizes roughness as measured by Sobolev semi-norm 
guarantees a maximum principle 

for the discrete Laplacian approximation 
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OTHER OPTIMAL TRIANGULATIONS 


Minimize max edge length - 0(n 2 ) Edelsbrunner, Tan 
Greedy Triangulation - Ofn 2 ) 

Minimum weight triangulation 
not known to be NP-complete 
not known to be solvable in polynomial time 
variant is NP-complete 
approximations used 


STEINER TRIANGULATION - RECENT RESULTS 

Chew (89) - Range: [30°, 1 20°] 

size optimal among all uniform meshes 
Baker, Grosse, Rafferty (88) - Range: [13°, 9(f] 
aspect ratio < 4.6 

Bern, Eppstein, Gilbert (90) - Range: [36°- 80 3 ] 
aspect ratio < 5 

Ruppert (93) - Range: [alpha, Pi-2 alpha] 

< aspect ratio < 

size optimal within a constant Caipha 
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HIGH ASPECT RATIO TRIANGULATIONS 


Delaunay triangulation can be unsuitable for high aspect 
ratio, body-conforming triangulations. 

Robust, efficient, global algorithms are in need. 

Computational geometers are not looking at this problem. 


SKEWED STRUCTURED GRID 



S 
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DELAUNAY REALIZABILITY 




DELAUNAY ANGLE CUT-OFF vs. ASPECT RATIO 
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ISSUES 

Generalize to a distance function which can 
vary throughout the plane. 

Avoid ambiguous cases. 
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CONVEX BODY PROJECTION AND CONVEX HULL 

Brown, 1979 
Edelsbrunner, 1987 

Project points from the plane to a paraboloid using parallel 
projection. 

Find the convex hull of the 3D point set (all points will be on the 
convex hull). 

The lower hull, projected back to the plane, will give the 
Delaunay triangulation of the point set in the plane. 

Notes: One convex body handles entire domain. 

Shifting the body to a new location gives the same result. 

CONVEX BODY PROJECTION AND CONVEX HULL 

Brown, 1979 
Edelsbrunner, 1987 
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STRETCHED TRIANGULATIONS 


Step la: Model simple stretching. 



STRETCHED TRIANGULATIONS 


Step 1b: Design convex surface which will produce 
desired stretched triangulation. 





STRETCHED TRIANGULATIONS 


Note: Body will not be "shift invariant 


paraboloid 


surface 
which gives 
stretched 
triangulation 
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Circumshapes derived from x 2 + y 4 , 6 — 0.02 



Triangulation derived from x 2 + y 4 
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Circumshapes derived from x 4 + y 4 , 8 = 0.02 



Triangulation derived from x 4 + y 4 
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Circumshapes derived from x 3 + y 3 , 8 = 0.09 



Triangulation derived from x 3 + y 





Circumshapes predicted from perspective projection, z pro j = -100, 6 = 0.05 



Triangulation derived from perspective projection, z pro j = —100 


308 


CONCLUSIONS 


Benefits of computational geometry - guarantees of 
grid quality 
efficient algorithms 

Many efficient triangulation algorithms are available, 
but high aspect ratio triangulations are not among them. 

Interdisciplinary cooperation will benefit grid generation 
and computational geometry. 


309 



